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ABSTRACT 

Q>^ ' We present an analysis of 3D spectra of Mrk 533, observed with the integral- field 

I spectrograph MPFS and using the Fabry Perot Interferometer (FPI) of the SAO RAS 

i 6-m telescope. We found emissions of gas from the active Sy 2 nucleus in the centre 

' and also from the HII regions in a spiral structure and a circumnuclear region. The 

, gas kinematics shows regular non-circular motions in the wide range of galactocentric 

distances from 500 pc up to 15 kpc. The maps of inward and outward radial motions 
of the ionized gas was constructed. We found that the narrow line region (NLR) is 
. composed of at least two (probably three) kinematically separated regions. We detect 

Qh' a stratification in the NLR of Mrk 533 with the outflow velocity ranging from 20-50 

' I kms~^ to 600-700 kms~^, respectively, on the radial distances of ^ 2.5 and ~ 1.5 kpc. 

^ . The maximal outflow velocity comes from the nucleus and corresponds to the position 

of the observed radio structure, which is assumed to be created in an approaching jet. 
^ [ We suggest that these ionized gas outflows are triggered by the radio jet intrusion in 

an ambient medium. 

. . Key words: galaxies: kinematics and dynamics - galaxies: Seyfert - galaxies: jets - 

' galaxies: individual: Mrk 533. 



1 INTRODUCTION 

Theoretical studies have shown that the interactions of 
galaxies can bring gas from the disc towards the nuclear 
regions, an d that they can also produce a burst of star 
formation l| Perez et al.l l2006h . Therefore, the detailed in- 
vestigation of the emission of gas kinematics in nearby 
Seyfert galaxies is very important. The panoramic (3D) 
spectroscopy is a powerful tool for this goal. Here we present 
a spectroscopic study of Mrk 533, an interacting Sy 2 
galaxy. Mrk 533 (also denoted as NGC 7674, Arp 182, 
UGC 12608, H96a) is the brightest member of th e Hick- 
son 96 compact group of galaxies (|Hicksonl Il982i ). It is 
the only spiral (Sbc pec) galaxy among other tidally in- 
terac ting galaxies in the group (jVerdes-Montenegro et ahl 
Il997h . The kinematics of Mrk 533 was observed twice 
with s canning Fabry - Perot systems l|Afanasiev fc KostiukI 
1 19961 : lAmram et al.l |2003[) and once using th e lon g 
slit-spectra technique ( Verdes- Montenegro et al.l Il997h . 

* E-mail: alexiya@sao.ru 



lAfanasiev fc KostiukI (|l996l ) mapped the velocity field in 
the Ha line observed with the 6-m telescope using FPI. 
The spectral resolution was S\ ~ 1.0 A and the free (in- 
terfringe) spectral range was AA « 13 A. They also used an 
integral-field data and reported a complex structure of emis- 
sion lines in the central region within r < 5 arcsec, where 
the Hq has a broad and narrow component. Moreover, the 
narrow component could be also divided into two Gaussians, 
where one of them corresponds to the normal rotation and 
secon d seems to be connect e d wit h the nuclear radio struc- 
ture. lAfanasiev fc KostiukI (|l996l ) found non-circular mo- 
tions that are proba bly caused by the circumnuclear bar. 
lAmram et al.l (|2003h also observed Mrk 533 and its envi- 
ronment with 3. 6-m ESO Telescope using the FPI system 
{SX « 0.7 A and AA « 8.5 A). They found a significant dis- 
agreement in shape of radial velocity distribution between 
receding and approaching sides of the galactic discs at dis- 
tances r > 20 arcsec. A possible shift of the rotation centre 
around 2 arcsec with respect to the optical centre was also 
noted. Both features are probably caused by interactions 
with companions, especially with HCG96c. 
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Mrk 533 is a Sy 2 galaxy where the broad Ha and 
H/j h nes are observed in polarized flux (jMiller fc GoodrichI 
Il990l ) indicating the presence of a hidden Broad Line Re- 
gion (|Tranl [l995l ) . In addition, the Narrow Line Region of 
Mrk 533 is extended, as s hown in the [C>III]AA 4 959,5007 im- 
age ob served with HST l|Schmitt et al.ll2003f ). rUnger et all 
1 19881 ) found that radio maps reveal the presence of a triple 
radio source with a total angular extent of about 0.7 arc- 
sec. This provides evidence that the radio emission is pow- 
ered by ejection. In the plane of the sky, the ejection axis is 
foun d to be roughly pe rpendicular to the galactic rotation 
axis (jUnger et al.|[l98a ). Recently, VLBI continuum and HI 
absorption observations of the central part with 100 mas 
resolution, showed six continuum structures extending over 
1.4 arcsec (742 pc), with a total flux dens ity of 138 mjy 
iMomiian et al.ll2003h . lMomiian et al.l l|200^ also suggested 
that the overall S-shaped radio pattern could be the result 
of a interstellar medium diverting the out-coming jets from 
the central AGN. In the same time, they could not role out 
the possibility of a black hole merger that could result in a 
similar structural pattern. 

The galaxy has noticeable blue-shifted wings on almost 
all of the emission li nes (except [SII] A A6716 ,6730A) as orig- 
inally suggested by lAfanasiev et al.l (Il980l) and later con - 



1981 ). 



firmed by other auth o rs (e.g . Shuder fc Osterbrockl 

Laurikainen fc MolesI ([1988'), 'P e Robertis fc Shawl ., . ., 

Voilleux (199^), etc.). Moreover, the observed UV doublet 
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O VIAA1032,1038 line profile shows a blue wing that ind i- 
cates an outflow in the U V emission gas (|Shastri et al.|[200^ . 
Also, the absorption lines at -300 kms ~^ and -800 kms~ ^ 
have been observed in the O VI doublet l|Shastri et al.ll2006h . 
The absorption in the O VI and a blue asymmetry in the 
optical lines might be caused by an accelerated wind from 
the central source. 

Here we present an analysis of Mrk 533 spectra observed 
with 6-m telescope of the SAO RAS using the integral-field 
spectrograph MPFS and Fabry-Perot Interferometer (FPI). 
The aims of the paper are: (i) to investigate the emission line 
structures of Mrk 533, focusing on the central part in order 
to map gaseous outflows indicated in the previous radio and 
UV observations; (ii) to study the gas kinematics in the cir- 
cumnuclear region. Also, here we will discuss the kinematics 
of a gaseous disc as a whole which has a non- circular com- 
ponent of motion in order to find the link between the host 
galaxy and central part kinematics. 

In this paper we adopt the value of 116 Mpc given by 
iMomiian et al.l ( |2003l ). to be the Mrk 533 distance, using the 
scale 1 arcsec=563 pc. 

The paper is structured as follows: In Section [2] we 
present our observations, in Section [3] we analyze the gas 
kinematics in Mrk 533, in Section |4] we investigate in more 
detail the complex NLR in this galaxy, in Section [5] we dis- 
cuss the possible link between the radio jet and optical one, 
and finally in Section [6] we give our conclusions. 




1 

-40 


-20 

-Aa (arcsec) 


1 

20 


1- 

40 


-200 


Velocities (km s"') 
-100 100 


200 


300 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




- , 1 , 
-40 


-20 20 
-Aa (orcsec) 


1- 

40 


20 40 


Velocity dispersion (km s"') 
60 80 100 120 140 


160 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




-20 

— Aa (orcsec) 



2 OBSERVATIONS AND DATA REDUCTION 

The panoramic (or 3D) spectroscopy provides spectra of ex- 
tended object simultaneously for each spatial element in a 



Figure 1. Large-scale emission-line intensity distribution in log- 
scale (a), ionized-gas line-of-sight velocity field minus the systemic 
velocity (b) and velocity dispersion (c) for Mrk 533 obtained from 
FPI data for the Ha line. The cross corresponds to the continuum 
centre. The thick contour on (b) marks the zero velocity. 
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Table 1. The spectroscopic observations of Mrk 533. 



Date 


T 


Sp. Range 


Sp. Res. 


Seeing 


Instrument 




(sec) 


(A) 


(A) 


(arcsec) 




1998 Aug. 30 


24 X 180 




3.0 


1.5 


FPI 


2002 Aug. 13 


3600 


4500 - 7100 


7.5 


1.5 


MPFS 


2005 Sep. 28 


6000 


4350 - 5900 


4.2 


1.3 


MPFS 



two-dimensional field of view (FOV). We used two differ- 
ent types of 3D spectrographs: the integral-field spectro- 
graph MPFS and scanning FPI. The first device provides 
a detailed spectral information in the spectral range about 
several thousand A what is important for the goals of the 
spectrophotometry and for analyzing the broad line profiles. 
However, the FOV of MPFS is small (less than 16 arcsec), 
it can cover only a region around the active nucleus of Mrk 
533. In contrast, the scanning FPI had a very large FOV (a 
few arcmin) with a very narrow spectral region containing 
only a spectral line. Using this device we are able to investi- 
gate the long-scale kinematics of the ionized gas in the whole 
stellar disc of Mrk 533, but only in the narrow component 
of the Ha emission line. 

2.1 Integral-Field Spectrograph MPFS 

Mrk 533 was observed in August 2002 and September 2005 
using the MultiPupil Fiber Spectrograph (MPFS) , with the 
integral-field u nit mounted at the primary focus of the 
6-m telescope (|Afanasiev. Dodonov fc MoiseevI l200ll ). The 
MPFS takes simultaneous spectra of 240 spatial elements 
(constructed in the shape of square lenses) that form on 
the sky an array of 16 x 15 elements (in 2002), or of 256 
spatial elements (16 x 16 in 2005). The angular size was 
1 arcsec per element in 2002, and 0.75 arcsec per element 
in 2005. The detectors were a CCD TK1024 (1024 x 1024 
px) in 2002 and EEV42-40 (2048 x 2048 px) in 2005. A de- 
scription of the MPFS is available at SAO RAS web page 
\http:/ /www, sao.ru/hq/lsfvo/devices. html\ 

During the first run, in August 2002, for spectrophotom- 
etry purposes, the galaxy was observed at a low resolution, 
while during the second run it was subject to a kinematics 
investigation of the NLR, and was observed at a higher res- 
olution. The log of MPFS observations is given in Table [T] 

The data were reduced using the software developed 
at the SAO RAS by V.L. Afanasiev and A.V. Moiseev 
and running in the IDL environment. The primary reduc- 
tion included bias subtraction, flat-fielding, cosmic-ray hits 
removal, extraction of individual spectra from the CCD 
frames, and their wavelength calibration using a spectrum of 
a He-Ne-Ar lamp. Subsequently, we subtracted the night-sky 
spectrum from the galaxy. The spectra of the spectropho- 
tometry standard stars were used to convert counts into ab- 
solute fluxes. 

2.2 Scanning Fabry-Perot Interferometer 

The galaxy was observed in August 1998 with the scan- 
ning FPI. The Queensgate interferometer ET-50 was used 
in the 235th interference order (for the Ha line). The free 
spectral range between neighboring orders (interfringe) was 



about 30 A. For monochromatization a narrow-band filter 
{FWHM — 3lA) centred on the spectral region contain- 
ing the redshifted Ha was used. During the observations 
we successively took 24 interferometric images of the object 
with different gaps between the FPI plates. The spectral 
channels were of SX « 1.2 A width, (~ 56 kms"^), the 
spectral resolution (the width of instrumental contours) was 
FWHM ^ 3.0 A (~ 140 kms"^). The detector was a CCD 
TK1024 (1024 x 1024 pixels) operated by 2 x 2 binning for 
reading-out time economy. Therefore each spectral channel 
had a 512 x 512 pixels format with the scale 0.68 arcsec 
per pixel. The log of the FPI observations is presented in 
Table [B 

To reduce t he interferomet ric observations we used a 
special software (|Moiseevll2003 ). running in the IDL envi- 
ronment. After the primary reduction (bias, flat-field, cos- 
mic hits) , we removed the night-sky spectrum, converted the 
data to a wavelength scale, and prepared them as a 'data 
cube'. The spatial resolution after data reduction smoothing 
was about 2.5 arcsec. The velocity fields of the ionized gas, 
and images in the Ha emission line were mapped using the 
Gaussian fitting of the emission line profiles. Moreover, we 
created an image of the galaxy in the continuum close to the 
emission line. 



3 THE KINEMATICS OF MRK 533 
3.1 Analysis of FPI observations 

Both previous observations were performed us ing the FPI 
(j Afanasiev fc Kostiuklll996l : lAmram et al.|[2003h with a high 
spectral resolution SX = 0.7—1.0 A but within a small spec- 
tral range (AA = 8.5 - 13 A). We performed new observa- 
tions with the 6-m telescope with a relatively low spectral 
resolution [SX ~ 3.0 A), but with a larger AA « 30 A. 
It is important for the study of the emission line profile in 
the nucleus, where FWHM of lines reaches 8 — 10 A. Un- 
fortunately, this spectral range is insufficient for a detailed 
analysis of emission line profiles, with FWHM larger than 
30 A (see Section gj. Therefore, with FPI we have consid- 
ered only kinematics of the narrow component of the Ha. 
For investigation, the broad component and the blue-wing 
asymmetry of the [OIII] lines we have used MPFS data only. 

The Ha monochromatic image, obtained with FPI, is 
shown in Fig.[T^. As one can see in Fig.[T^, besides the strong 
central part in the Ha, there are structures (which follow the 
bar and spiral arms) which show a significant Ha emission. 
The ionized gas velocity field and velocity dispersion map 
(corrected on the instrumental contour) are shown in Figs 

As can be seen in Fig. [TJj, the rotation axis is SW- 
NE oriented. It is interesting to note that in the velocity 
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Figure 2. Results of the tilted-ring analysis of the H« velocity 
field: radial variations of the rotation velocities (top) and position 
angles (bottom). Full circles, stars and open diamonds denote the 
value calculated on the whole disc, SE (approaching) and NW 
(receding) sides respectively. Dotted line corresponds to PAo = 
306° adopted in this paper (see for details Section 13.21 . 

map two features with approaching velocities (with respect 
to the systematic one) are present. These features will be 
discussed in more details in Subsection 14.41 The velocity 
dispersion map (Fig. [!};) shows the arms near the centre 
and, as expected, the maximum velocity dispersion in the 
central part. 

Our maps seem to b e slightly deeper than the ones given 
by lAmram et al] (j2003l ). For instance, we detected the ve- 
locities in several emission clumps along tidal tail located 
at the North of the galaxy. The field-of-view contains also a 
nearby satellite HCG96c located to the NE from the Mrk 
533. We measured the systemic velocity of this compan- 
ion, and obtained 8850 kms~^ (z =0.0295). This is in good 
agreem ent with the value given bv lVerdes-Montenegro et al.l 
I^1997^■ and diffe r ent a round ~ 200 kms~^ from the results 
of lAmram et all (120031 ). 

3.2 The model of pure circular rotation 

The velocity fiel d was ana l yzed using the so called 'tilted- 
rings' method bv lBegemanI i|l989l fl. We were looking for the 
parameters of the galactic gaseous disc orientation, i.e. a 
line of node position angle (major axis) PAq and the galac- 
tic plane inclination iq. The velocity field was fitted using 
a pure circular rotation model by fixing the centre posi- 
tion and systemic velocity. We excluded from the analysis 
the central region (r < 7 arcsec) where strong non-circular 
motions are present (see Section |4]) and the external points 

^ See also iMoiseev. Valdes fc ChavushvanI l|2004) and references 
therein for more details. 



(r > 25 arcsec) where the vel ocity field is significan tly asym- 
metrical, as already noted by lAmram et all (|2003l ). The fol- 
lowing values were found from a x^-minimization routine of 
residuals between the observed and modeled velocity fields: 
io = (33 ± 5)°, PAo = (306 ±4)°. We assumed that the 
position of the kinematic centre corresponds to the veloc- 
ity field centre symmetry. The obtained systemic velocity 
is Vsys = 8672 ± 3 kms"^ (z=0.0289). The kinematic cen- 
tre of the galaxy is shifted for f pixel (or 0.7 arcsec) to the 
South from the continuum centre. This centre offset (lo psid- 
cdne ss) was already ment i oned b y lAmram et al.l l|2003h and 
Vcrd es-Montenegro et all (| 19971 ). however different authors 
found difi'erent positions of the centre offset. The observed 
shift of the kinematic centre might be due to the interaction 
with the companion as well as with matter that has a non- 
circular motion caused by the bar's gravitational potential. 

The orientatio n parameters are in accor d ance with the 
ones measured by IVerdes-Montenegro et al.l (|l997l ). PA = 
121° (i.e. 301°), i = 31°, while these values are signi ficantly 
different to the ones reported bv lAmram et al.l(|2003l ). PA = 
132 ± 5° and i = 50 ±5°. 

3.3 Detailed analysis of the Ha velocity field 

In order to study the nature of the deviations from the nor- 
mal rotations, the velocity field was divided into elliptical 
rings with a 1 arcsec width aligned to the accepted PAq 
and io. Using the minimization routine for deviations of 
the observational points from the model, we calculated, for 
each radius, the position angles of the major axis PAdynir) 
(dynamical position angle) and the mean of rotational ve- 
locity Vrot{r). The results of the best-fitting are shown in 
Fig. [21 where results for approaching (asterisk in Fig. ^ and 
receding sides (diamonds in Fig. [2)) are shown. Th e se are 
calculated separately. As was noted in lAmram et all (|2003| ) 
earlier, we obtained a different kinematics between these two 
sides. 

The mean ('in the whole disc') rotation curve shows a 
very unusual shape at distances larger than r = 20 arcsec. 
The maximum is at r ~ 26 arcsec, after that the rotational 
velocity decreases. As it is shown in Fig. (2] this maximum in 
the mean rotational curve is caused by increase of the rota- 
tion velocity (after r > 20 arcsec) in the approaching (SE) 
side of the disc, while Vrot{r) of the receding side is slightly 
decreasing in this region. Also, a small Vrot{r) peak in the 
approaching side can be noticed between 10 and 14 arcsec. 
These features cannot be caused by specific mass distribu- 
tion, but rather by an interaction with other members of the 
Hickson group, especially at the maximum r > 20 arcsec. 

There are two main possibilities, both triggered by a 
close interaction: (i) non-c ircular gas motion s (first of all in- 
wards gas streaming, see llono et al] (|2004l ) and references 
therein) and/or (ii) quasi-circular gas motions with the ori- 
entation parameters which are differe nt from t h e inne r galac- 
tic disc ones (warped outer disc, see iBinnevI l| 19921 ) for re- 
view) . 

In case (i) we can estimate the peculiar velocities via 
the map of residual velocities (the differences between the 
observed line-of-sight velocity field and the circular model, 
see Section [3. 4p . The 'blue' excess of the observed velocities 
(—(30 — 60) kms~^) appear near the major axis at r > 25 
arcsec on both sides of the centre if we compare the observed 
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Figure 3. Ionized gas kinematics in the central region: (a) residual velocities after subtraction of the model of pure circular rotation 
from the observed velocity field, Ho isophote are overlapped, the box corresponds to MPFS field from Fig. llOb : (b) the map of radial 
motions: the inflow is denote with blue and the outflow is red. The galactic bar and spiral arms are shown. The central region mapped 
with MPFS (Figs. I5I9I I is shown as a small box in the centre. 
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Figure 4. Difference in the blue part of the [OIII] lines across the central (nucleus) part of Mrk 533 (a). Decomposition of the [OIII] lines 
with two (b) and three (c) Gaussian components. Dots represent the observation and the solid line represent the best-fitting. Gaussians 
are shown below. 



data with the model calculated for the whole disc (the full 
circles in Fig. [2}. On the SE side of the disc the peculiar 
velocities located in the outer part of the spiral arm are 
directed towards the companion HCG96c. An asymmetric 
shape of this part of the sp iral structure is due to the inter - 
action with the satellites l|Verdes-Montenegro et al1ll997^ . 
The line-of-sight projection of the radial motions is zero to- 
wards the major axis; using equation ([1} from Section |33] we 
can obtain the residual velocities in the SE spiral arm. That 
correspond to the vertical velocities of up to 70 kms~^ or 



to the azimuthal streaming in the galactic plane of 80 — 120 
kms~^ (in the direction of the rotation). These strong non- 
circular velocities are similar to the ones observed in the 
grand design spiral interaction galaxies, as e.g. in M51, wit h 
velocities ±(70-130) kms~^ reported bv lVogel et aD (jl993h . 

On the other hand, the situation on the NW side of 
the disc is different. The residual velocities in the outer re- 
gions of the NW disc side have similar amplitude as on the 
SE side, but it is only in the interarm region and it has no 
connection with the morphological structure of the disc. At 



6 A. A. Smirnova et al. 



the same time the deviations from the circular model are 
around zero in the NW spiral arm and in the tidal tail at 
the North. Moreover, if we calculate the model of rotation 
only for the NW side of the disc (the diamonds in Fig. [2]) 
the peculiar velocities in the interarm region decrease sig- 
nificantly. Therefore, the non-circular velocities here seem 
to be linked to the global distinctions in the orientation pa- 
rameters of the gaseous disc at a different values of r (i.e. a 
warped disc). 

If there is a case of warped disc, one can estimate the 
warping angle of the disc assuming a constant rotation ve- 
locity along the radius. We found that for r — 35 arcsec the 
inclinations are i' = (20 ±4)° and i' = (48 ± 5)° for SE and 
NW sides of the disc, respectively. Therefore, if the disc is 
warped between r — 20 arcsec and 35 arcsec is moderate, 
\i' - io\ = (13 ± 6)° for SE and (15 ± 7)° for NW sides. As 
we already mentioned above, it is possible for both mecha- 
nisms to contribute to the observed kinematics: a significant 
non-circular motions along the SE arm and warped disc in 
the NW side of the galaxy. 

The radial variations in PAdyn around the line-of-nodes 
is 10 — 20° in amplitude (see Fig. 01 ■ It indicates non- 
circular gas motions. Exactly, at r < 7 arcsec the dynam- 
ical axis turns in the opposite side to the major axis of 
the near-infrared isophote of the circu mnuclear bar wit h 
PAphotom^^try « 280° - 290° (scc Fig. [T]in |H^^illL|[l993) . 
This type of PAdyn variation corresponds to the radial 
streaming motion along the b ar (see discussion and ref- 
erences in iMoiseev et al.l 120041 ') . We note here that these 
motions should be bi-symmetrical, since measurements of 
PAdyn show a good agreement between the SE and NW 
sides of the disc. Comparise to the bar's region, the PAdyn 
deviations at large distances (r = 13 — 22 arcsec) are signif- 
icantly different at the SE and NW sides of the disc. The 
possible explanation to this mismatch, is an influence of the 
spiral arms on gaseous dis c kinematics. The shape of t he sp i- 
ral pattern is asvmmetric IVerdes-Montenegro et al.l ([1993), 
hence the Northern and Southern arms cause different am- 
plitudes (and characters) of velocity perturbations, as it is 
shown in Fig. [5] 

3.4 The map of non-circular motions 

For a precise analysis of the gas kinematics in the circumnu- 
clear region, we study a spatial distribution of possible non- 
circular motions. With this goal we made a map of residuals 
between the observed velocity field in the Ha and the model 
of the pure circular rotation. The circular model was ob- 
tained using the 'tilted-rings' approximation (see previous 
section) with fixed angles PA — PAq and i = io. Fig. [5^ 
shows a 10 X 10 kpc region of the residual velocity field. This 
figure reveals a tight correlation between large amplitude 
of the residuals and spiral arms traced in the Ha-isophote. 
Moreover, large non-circular velocities are detected in the 
centre (r = 1 — 5 arcsec) of the galaxy along its minor axis: 
Vres = +{20 - 25) kms"^ to NE and Vres = -(40 - 45) 
kms~^ to SW from the nucleus. Before interpreting these 
features, let us recall a formula for the projection of veloc- 
ity vector in the galactic disc after the subtraction of the 
circular rotation: 

Vrea = Vr slu tfi siu i + V^p COS siu 1 + Vz COS i, (1) 




Figure 5. The outflow velocity map (left) and the field of ve- 
locity dispersion (right) of the NLRl, where assumptions of two 
emission regions were applied. The zero-point cross marks to the 
centre in the continuum images. 



where {Vr,Vip,Vz} are radial, azimuthal and vertical compo- 
nent of the velocity vector respectively; ip is the azimuthal 
angle measured from the projected major axis in the plane 
of the galaxy. Assuming that the motion is only in a galactic 
plane, than according the first approximation, we have that 
Vz ~ 0. Therefore, in the points which are close to the minor 
axis {ip « 90,270°) we observe only a projection of radial 
motions: Vres ~ Vr sini. The direction of these motions (in- 
flow or outflow) is determined by the orientation of the disc 
relative to the o bserver. For Mrk 53 3 the SW side is closer 
to the observer iKinnev et al.ll2000l ). that is in accordance 
with the assumption of trailing spiral arms. In this case, 
large residual velocities along the minor axis correspond to 
the outflow motion 40 — 90 kms~^ in an amplitudsiJ. 

Using this formalism we created a map of inflow/outflow 
radial streaming motions (Fig. [3]d). On this map we avoid 
points I V' l< 45° that form with the major axis. Here the 
projection of the azimuthal component of residual velocities 
can be larger than radial ones (for this points should be 
I sin </3 I < I cos ifi I because in spiral galaxies Vr and V^ have 
same amplitude, see direct measurements by Fridman et al. 
( 2001a . bi). Also we show in Fig. [Sja the position of the bar 
(with orientation and axis ra tio is in accordan ce with the 
NIR phot ometry reported bv | Hun t et al."l997l) and spiral 
arms fsee IVerdes-Montenegro et al . (1997), their Fig. [Gja). 
The map of radial motions shows complex structures: the 
outflow in the centre (r < 5 arcsec), the gas inflow along the 
Northern arm at r = 5 — 15 arcsec and an outflow again at 
r > 15 arcsec (8.5 kpc) along the Southern spiral arm. The 
difference in radial motions for the different arms can be a 
result of their asymmetrical shapes caused by an interac- 
tion with other members of the compact group of galaxies. 
Also the possible position of a co-rotation at r ~ 15 arcsec 
may produce these changes of radial motions in the direction 
along the spiral arm. 



^ This speed should decrease if the direction of the outflow lies 
out of the galactic plane, i.e Vz 7^ 
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Figure 6. The outflow velocity map (left) and the field of ve- 
locity dispersion (right) of the NLRl, where assumption of three 
emission regions was applied. 



4 THE NLR KINEMATICS 

In order to investigate the NLR kinematics, in this section 
we present the analysis of the [OIII] emission Une profiles 
obtained from MPFS data. The shapes of the narrow spec- 
tral lines have a blue asymmetry only in the nuclear part 
(see Fig. |4]) of the galaxy. By inspecting the line profiles, the 
asymmetry can be seen in the Ha, H/3 and [OIII] lines, and 
the shapes of the blue part of the lines differ from the various 
parts of the nucleus emitting region (Fig. [Jjl . Here we have 
high spectral resolution spectra only for the H/3-|-[0III] lines, 
observed in 2005, but the H/3 line is significantly weaker than 
the [OIII] lines and asymmetry is not enough pronounced. 
Consequently we analyze in more details only the spectra of 
the [OIII] lines. 

It is important to note that we could spatially resolve 
the region with a blue-wing asymmetry of line profiles. In 
order to check this, we created a map showing the distri- 
bution of brightness of the fiux in the broad wing of the 
[OIII] lines. Indeed, our analysis showed that the FWHM 
of the spatial distribution was about 2-2.4 arcsec that is 
significantly larger than the size of the image of the star 
(seeing was 1.3 arcsec) observed with MPFS directly before 
the galaxy. 



4.1 Line profile analysis 

To analyze the shape of the [OIII] lines one can fit the line 
profiles assuming different functions (Lorentz, Voigt, Gaus- 
sian profiles). But taking into account that the line profiles 
in the emitting region is caused by a motion of the emitters, 
here we fitted each line with the sum of Gaussian compo- 
nents. We applied a minimalization routine in order to 
obtain the best-fitting parameters. We also assumed that 
the narrow emission lines can be represented by a sum of 
Gaussian components. In the fitting procedure, we looked 
for the minimum number of Gaussian components needed 
to fit the lines. To limit the number of free parameters i n 
the fit we set some a prior constrains (jPopovic et al.ll2004h : 
(i) Two Gaussian (e.g. broad and narrow components 
shown in Fig. [5]) representing the [OIII]AA4959,5007 lines 
are fixed at the same red-shift and the Gaussian widths are 



Figure 7. The outflow velocity map (left) and the field of velocity 
dispersion (right) of the NLRla. 
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Figure 8. The field of velocity dispersion of the NLR2: 1. The 
case of fitting spectra with two Gaussians (left); 2. For the case 
of fitting spectra with three Gaussians (right). 

set proportional to their wavelengths, 

W4959 _ W5007 
4959 5007 

(ii) We assumed that the i ntensity ratio of the t wo 
[OIII] AA4959, 5007 lines is «1:3 (|Dimitriievic et al.ll2007l ). 

Here we will use the velocity dispersion to present the 
velocity field (a = {W/\/2) ■ (c/A), where W is the Gaussian 
width). Outfiow velocities are considered to be K- = c x 

where A A = Ac — Af, is a blue shift, measured between 
the centre of the narrowest component that is assumed to 
correspond to the transition wavelength Ac and those which 
are blue shifted At (Fig. H]) 

As mentioned above (see in Fig. |4} the [OIII] lines show 
a blue asymmetry, and can be properly fitted with two Gaus- 
sian components (Fig. HJj). Moreover, a slightly better fit 
(concerning x^) can be obtained by fitting with three Gaus- 
sians (Figll]:). Therefore, here we will consider the complex 
[OIII] line region as composed from two or three separated 
emission regions. When we assume the emission of two re- 
gions, we denote those regions as the NLRl (blue shifted 
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Figure 9. Flux ratio between NLRl and NLR2 emission com- 
ponents where two (left) and three (right) emission regions are 
assumed. 



component) and NLR2 (non-shifted component). In the case 
where we assume that the emission is coming from three re- 
gions, we additionally consider a NLRla existence (also the 
blue-shifted component) region that corresponds to a com- 
ponent placed between the central and broad blue shifted 
one (Fig. 3};) . We started to fit each spectra obtained from 
spatial MPFS elements, but out of the nucleus (around 8x9 
elements) it is not possible to see the asymmetry, so the 
method of fitting analysis has been applied only in the small 
frame elements that corresponds to a small region on the sky 
around 6x7 arcsec. 

The velocity maps that we obtained from our analysis 
are presented in Fig. 1 5181 Also we found a map of flux ratios 
of NLRl and NLR2, i.e. I^ii-Si (Fig.©. 

4.2 Assumption of two emitting regions 

Assuming the presence of two kinematically separated re- 
gions, we fltted the [OIII] lines with two Gaussians (Fig. 
|4j3). As one can see in Fig. |4]3 the line shapes can be well 
fitted with two Gaussians, but a small shoulder in the blue 
part can not be fitted properly with this assumption. The 
map of the outfiow in the NLRl is presented in Fig. [S]left. 
As can be seen in this figure, the maximum outflow velocity 
(~ —550 kms~^) is located to the east from the contin- 
uum centre (presented as a cross in Fig. [5} ■ A gradient in 
outfiow velocities from the east to the west is present. The 
obtained difference amongst velocities (from the east to the 
west) is ~ 300 kms~^. This supports the idea about the ex- 
istence of an approach ing jet from the active galactic nucleus 
(jMomiian et al.ll2003l ). The velocity field of the outflow has 
a south-north elongated structure (Fig. [5^), which might be 
due to the projection of the jet. On the other hand the veloc- 
ity dispersion field has a relatively small difference in NLRl 
(~ 560-670 kms~ ), showing a relatively flat velocity fleld. 
The intensity ratios in the central part of NLR1/NLR2 are 
present in Fig. (left). The maximal intensity of the NLRl 
is around the optical centre. 

The NLR2 velocity dispersion map shows also stratifica- 
tion, where velocity dispersion has values ranging from 150 
kms~^ to 170 kms~^, having a gradient from the west to 
the east (Fig.[8]left). This gradient in velocity dispersion has 



a possible connection with the radio jet which is also elon- 
gated in the same direction (see Section [5]). However this 
should be taken with caution, since the differences in veloc- 
ity dispersion are too small (around 20 kms~^, that is very 
close to the range of the error bars of the fitting method). 

4.3 Assumption of three emitting regions 

As was mentioned above, a slightly better fit can be ob- 
tained when we perform a three Gaussian fitting procedure. 
Moreover, there is an anticorrelation between the outflow 
velocity map of the NLRl and dispersion velocity map of 
NLR2 (see Figs. [5] and [Si that may be caused by fltting the 
observed profile with an insufficient number of components 
(only two components). Therefore, we assume that in the 
nucleus of Mrk 533 three, kinematically separated, emission 
regions are present on the line-of-sight. The velocity maps 
of the NLRl and NLRla are presented in Figs. |6] and [T] 
As one can see in Fig. [6] the outflow and velocity dispersion 
in NLRl are slightly higher then in the case of two region 
approximation. The maximum outflow velocities (^ —700 
kms~^) are located to the east of the optical centre, hav- 
ing a south to north elongated structure. In comparison to 
the two region approximation, the three region approxima- 
tion shows that the outflow structure is more complex then 
the one obtained from two Gaussian fit. Also, the map of 
velocity dispersion shows several regions around the optical 
centre with higher velocities (Fig. [6l right). But this should 
be taken with caution, since the difference in velocity disper- 
sion is small (~ 100 kms~^) and the velocity dispersion field 
of the nucleus of Mrk 533 tends to have a flat distribution, 
as it was the case in the two Gaussian assumption. 

The NLRla region shows a complex structure around 
the optical centre, with maximum outflow velocity around - 
250 to -300 kms-\and minimum velocity dispersion around 
140 - 320 kms"\ see Fig.[71 Also, the ratio between NLRl 
and NLR2 fluxes (Fig. ^ indicates a clumpy NLR. 

In this case the NLR2 velocity dispersion region is simi- 
lar to the one in two Gaussian case, where velocity dispersion 
has values ranging from 70 kms~^ to 160 kms~^ (Fig[8]). As 
one can see there are some structures which might indicate 
also a clumpy NLR2. 

4.4 The complex NLR: Two or three emitting 
regions? 

Based on our analysis the [QUI] line shapes across the cen- 
tral part of Mrk 533, it can easily be noted that there is an 
outfiow. The outflow also can have a complex structure in 
velocity and spatial domain. Although the [OIII] line pro- 
flies can be well fitted with the assumption of two emitting 
region (where one is characterized by the presence of an out- 
flow) , the better flt of the line shapes can be reached with 
the three emitting region assumption. Moreover, the anti- 
correlation between the NLRl outflow velocity map and the 
NLR2 velocity dispersion map in the two emission region 
indicates that there are probably two regions with outflows. 
Here the outfiow velocities range from -500 to -700 kms~^ 
in one region and range from -200 to -400 kms~^ in the 
another. The obtained velocities are in accordance with the 
outfiows registered from the O VI absorption lines, where the 
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measured velocity o f the outflows are -800 and -300 kms^^ 
IShastri et al.ll2006h . It seems that these regions in outflow 
also absorb in the UV radiation. On the other hand, our re- 
search of the circumnuclear region (see Section [S} shows that 
near the nucleus r < 5 arcsec there is an outflow, but with a 
smaller outflow velocity (from -40 to -90 kms~^). This type 
of outflow was observed only in the narrow component of the 
Ha profile when the circular rotation velocities were taken 
into account. 



5 CONNECTION BETWEEN THE RADIO, UV 
AND OPTICAL OUTFLOW 

In this part we compare our results obtained in the optical 
range with the previous results which were obtained in the 
UV and radio band in order to connect the outflow(s) seen 
in the NLRl with the radio jets. An indication of a jet- 
cloud interaction is the ratio of the emission lines (which can 
originate from shocks), therefore we first explore the nature 
of the gas ionization source in the nuclear and circumnuclear 
region as well as in the spiral structure. 

5.1 The source of ionization 

In order to identify the nature of the gas ionization source in 
Mrk 533 we constructed diagnostic diagrams using the emis- 
sion line intensity ratios, i. e. in Fig.lTOlwe present [ QUI] /H/3 
vs. [NII]/Hq. According to lVeilleux fc Osterbrockl l|l987l ) we 
separated the regions corresponding to the ionization by 
AGN, young OB-stars (HII regions) and shock waves (LIN- 
ERS) into presented diagram. As a result of short analysis of 
the spectra, we were able to separate several emitting regions 
(see Fig. 1101 right). Mainly there are three emission regions 
in the stellar disc plus an emitting region that corresponds 
to the nucleus. The analysis of the flux ratio of [OIII]/H/3 
vs. [Nil] /Ha lines shows that spectra belonging to the nu- 
cleus (crosses) are located in the part where the ionization 
of AGN is expected. On the other hand, spectra from the 
spiral- like emission structure (regions 1, 2 and 3) are mainly 
located in the HII part of the diagram, i.e. the ionization 
from hot young stars is present. Note that some points from 
the region 2, that is closest to the AGN, are located also in 
the AGN part of the diagram. It means that non-thermal 
ionization predominates in some parts of the region 2. It is 
interesting that some points from region 3 are located near 
the HII/LINER borderline. This may indicate that besides 
the star bursts ionization, gas may be compressed and that 
in parts of this region, a shock ionization occurs. 

From the inspection of Fig.[TO]it is evident that different 
ionization mechanisms may dominate on different scales. We 
flnd that in nucleus part as well as in the regions (1 and 2) 
which are close to the nucleus the dominant ionization mech- 
anism is AGN-photoionization. However, as we move further 
away from the nucleus (~ kpc), jet-cloud interactions may 
become the main ionization mechanism. In order to discuss 
the possibility of the dominance of a jet-cloud interaction 
in some em ission regions of Mrk 533, let us recall results 
obtained bv lDopita fc Sutherland! (Il995l ). They investigated 
the line fluxes expected from fast, radiative shocks prop- 
agating through typical narrow-line AGN clouds. Here, the 
key parameters controlling the emission-line ratios are shock 



veloci ty and the magnetic parameter. iDopita fc Sutherland! 
(| 19951 ) found that inclusion of 'precursor' ionization caused 
by hard radiation generated in the wake of the post-shock, 
cooling plasma that diffuses outwards, is crucial for explain- 
ing differences found in Seyfert and LINER galaxies. It is 
interesting that faster shocks are required (~ 300 - 500 km 
s~^), that is well fltting our obtained outflow velocities in 
the NLRl and NLRla. We present in Fig. [10] with fuU cir- 
cles the predicted line ratio s for different shock velociti es 
('shock plus precursor', see iDopita fc Sutherland! l|l995h ). 
The observed [OIII]/H/3 vs. [Nil] /Ha ratios in the nucleus 
region are consistent with a 'shock plus precursor' model 
that has shock velocity exceeding 300 km s~^ at least, and 
a magnetic parameter typic al for the interstellar medium 
(jPopita fc Sutherland! I1995I ). Moreover, some points from 
the other three emission regions shown in Fig. llOb are also 
consistent with a 'shock plus precursor' model, but with 
smaller velocities ~ 200 km s~^. Note here that the flux ra- 
tios of [0III]/H/3 and [Nil] /Ha are constructed only using 
narrow lines observed with the MPFS. 

5.2 Outflovirs in the radio, optical and UV 
emission 

To reconstruct the picture obtained from the analysis of the 
NLR and circumnuclear kinematics we draw Fig. 1111 where 
the outflow velocities in all three components are given. The 
velocities are measured along the dashed line shown in Fig. 
1121 The circles in Fig. [TT] represent the measured velocity 
of the outflows from O VI absorption lines (-800 and -300 
kms~^, see lShastri et al.!l2006l ). As it can be seen in Fig. 1111 
there is a stratification in the line-of-sight velocities, that 
may also indicate a different distribution of these compo- 
nents, as well as, heaving in mind the different absorption 
densities of the UV. The spatial structure of the outfiow 
is slightly resolved: the high- velocities of {Voutfiow > 100 
kms~^) NLRl extends up to 1.5 kpc, and the low-velocity 
{Voutfiow = —(20 — 80) kms~^) emission regions extend up 
2.5 kpc (see Fig. Ilip . We remind that our spatial resolutions 
were about 0.7 kpc (1.3 arcsec) on MPFS and 1.4 kpc (2.5 
arcsec) with FPL Since the emission originates far away from 
the nucleus, one can consider the possibility that the mech- 
anism for the optical line emission may actually be a shock 
wave, created in a wind (or jet) from the AGN, that is also 
in accordance with the analysis of the source of ionization 
(see Section [3. 4|) . 

Note here, that with the MPFS and FPI we have de- 
tected essentially different components of the outflow. In the 
first case the multicomponent fltting of the emission lines 
profiles was carried out and the high-velocity component in 
the blue wing of the emission lines profiles was studied. In 
the second case, the narrow spectral range of the FPI did 
not allow us to analyze the shape of the spectral profile. 
However, higher (in comparison with MPFS) spectral reso- 
lution gives the most accurate measurements of the veloci- 
ties of the emission line core. Moreover, using a large field 
of view and better spatial sampling we can determine the 
differences between the observed velocities and the circular 
rotation model. Therefore, we can distinguish the peculiar 
velocities with amplitudes about several tens of kms~^. We 
illustrated the differences between the MPFS and the FPI 
measurement on the right side of Fig. [Tl]where components 
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Figure 10. The flux ratio of [OIII]/H/3 lines vs. [NII]/Ho (a) and the map of Mrk 533 with regions for w hich this ratios are found (b) . 
The Ha isophotes are overlapped on this map. The full circles in Fig. 4a represent the calculated ratio bv lDopita fc Sutherlaiidl lll995l) . 
where 'shock plus precursor' model were used. The calculated values are for the shock velocities of 200 km s~^, 300 km s~^, and 500 km 
s~^ (the values are given near the full circles). 



of the emission lines correspond to the given outflow veloc- 
ities shown schematically. 

To find any connec tion between the op t ical and the ra- 
dio ou t flow (reported in Unger et al.l () 19881 ) , iMomiian et al.l 



l|2003t ). lMiddelberg et a"l.ll2004l ) we draw a map with the op- 
tical and radio outflow structures (see Fig. I12p. It is indica- 
tive to see that the orientation of two jets from two points 
such as radio sources, are directed to the optical outflow 
structures. Also, as one can see in Fig. 1121 the maximal out- 
flow velocities correspond to the side where the most inten- 
sive part of the radio structure is located. Concerning the 
map shown in Fig. 1121 one can speculate that the optical 
outflow is following the radio one, i.e. that the jet emits in 
the radio domain, as it is close to the central engine (r < 1 
arcsec or 0.5 kpc), and farther from the centre of the UV ab- 
sorption/emission (see Fig. Ilip and optical emission occurs 
(r = 2 — 4 arcsec or 1 — 2 kpc). We outline that this opti- 
cal emission, partly triggered by the shock waves originates 
from the jet intrusion in to the gaseous ambient medium. 

The presence of high- velocity clouds, located close to 
the passage of the radio jet and centred near the centre of 
the nucleus, shows that the effects of jet-induced shocks are 
still important in the nuclear regions. Moreover, the veloc- 
ity dispersion of the [OIII] component in the NLRl is high: 
ranging from 570 to 720 km s~^. Increased linewidths are 
expected following turbulent motions of an interaction (e.g. 
by jet-cloud interaction), but not in passive photoioniza- 
tionby the AGN (see e.g. lClark et al.1 (|l997l ): iGandhi et al.l 
l|2006t )). On the other side in both approximation (two and 
three emission line region) we noted outflows with a rela- 
tively high velocity. Therefore, the assumption that the pro- 
jected radio jets (see Fig. I12|) have an influence on NLRl 
physics and kinematics seems to be real. 



6 CONCLUSION 

In this paper we analyzed the spectroscopic observations of 
Mrk 533 performed by the 6-m telescope of SAO RAS. We 
investigate the gas kinematics in the spiral and circumnu- 
clear structures. Based on this investigation we are able to 
conclude the following: 

i) We found the parameters of the galactic gaseous disc 
orientation: io = (33 ± 5)° and PAo = (306 ± 4)°. The 
parameters are in accordance with the ones obtained by 
IVcrdes-Montcnegro ct al. (1997), bu t they are significari tly 
different from the ones measured bv lAmram et al.l l|2003l ). 

ii) Our analysis of the large-scale ionized gas velocity 
field shows presence of significant non-circular motions. Also 
we indicate the outer warping gaseous disc, probably caused 
by an interaction with the companions. The disc warping 
was found to be around 13— 15° at distances r > 20 kpc (at 
least in the NW side of the disc). 

iii) We created a map of radial motions in the gaseous 
disc. It shows a very complex non-circular gas motion in 
Mrk 533, that may be triggered by an interaction. So, the 
gas inflow along the Northern arm (5" < r < 15") and 
streaming motions along the bar were found by analyzing of 
this map and as well as radial dependence of the dynamical 
axis PA. A co-rotation at r « 8.5 kpc may produce changing 
of radial motions direction along the Northern spiral arms 
at r > 15 arcsec. The circumnuclear outflow in r < 5 arcsec 
along the minor galactic axis was also detected. 

iv) The central emitting region of Mrk 533 is complex 
and can be divided into: the emission of the Sy 2 nucleus (in 
the very central part) and the star-forming regions which 
follow the spiral structure of the galaxy. The emission line 
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Figure 11. The outflow velocities in the central part of Mrk 533 
(within —4.5" < r < 5") obtained from Gaussian fits of the [OIII] 
lines (nuclear region) and from the analysis of circumnuclear kine- 
matics of Ha (full circles, crosses indicate the interval of scattering 
^outflow S'lid r). The circles represent t he outflow velocitie s from 
absorption components in O VI lines l lShastri et al.l|2006|) . The 
diagrams on the right show a component of the emission profile 
(marked as black Gaussian) corresponded to the different scales 

of Voutflow 

shapes in the region of the nucleus have a blue asymmetry 
this indicates an outflow from the centre. 

To investigate the velocity structures of the nucleus, 
we assumed that the NLR of Mrk 533 is composed of: (a) 
two and (b) three kinematically separated regions. From this 
analysis we can conclude that: 

- The [OIII]AA4959,5007 lines originate from at least 
two kinematically separated regions: NLRl contributes to 
the blue part of the spectral lines, where a gaseous outflow 
with velocities from -300 kms"'^ (5 arcsec far from the cen- 
tre) to -600 kms~^ around the nucleus (the velocity disper- 
sion is around 560 - 670 kms~^) exist; NLR2 has a smaller 
velocity dispersion (~ 150 kms~^). 

- The complex [OIII] lines may be produced in three 
kinematically separated regions: NLRl contributes to the 
blue part of the spectral lines, with a gaseous outflow with 
velocities ranging from -500 to -700 kms"^ the NLRla with 
outflow velocities from -350 to -100 kms~^ and NLR2 with- 
out a significant outflow. 

Finally, we can conclude that a stratification in the NLR 
of Mrk 533 exists, and that it has a range of outfiow veloc- 
ity from 20-50 kms~^ to 600-700 kms~^ on radial distances 
~ 2.5 and ~ 1.5 kpc accordingly. The structure of the opti- 
cal outfiow seems to be complex, but it seems that the reg- 
istered outfiow(s) in the optical corresponds to the one(s) 
registered in the UV. Moreover, it is likely that the outfiow 
obtained from the optical emission lines ([OIII] and Ha) fol- 
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Figure 12. The optical vs. radio outfiow: Gray colors shows the 
position of circumnuclear outflow region seen on Fig.lSb. The iso- 
lines show the outflow mapped using the [OIII] lines (Fig. |6] left) 
and feature in the centre is the map of the radio continuum given 
by l|Momiian et al.l l2003l l. The dashed line marks the direction of 
the X-axis on Fig. 1111 



lows the radio one and that the complex structure of the 
optical outfiow is generated by the radio jet. 
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